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The finite range adiabatic wave approximation provides a practical method to analyze (d,p) or 
(p,d) reactions, however until now the level of accuracy obtained in the description of the reaction 
dynamics has not been determined. In this work, we perform a systematic comparison between the 
finite range adiabatic wave approximation and the exact Faddeev method. We include studies of 
n Be(p,d) 10 Be(g.s.) at E v =5, 10 and 35 MeV; 12 C(d,p) 13 C(g.s.) at E d =7, 12 and 56 MeV and 
48 Ca(d,p) 49 Ca(g.s.) at Ed =19, 56 and 100 MeV. Results show that the two methods agree within 
« 5% for a range of beam energies {Ed f» 20 — 40 MeV) but differences increase significantly for very 
low energies and for the highest energies. Our tests show that ADWA agrees best with Faddeev 
when the angular momentum transfer is small Al = and when the neutron-nucleus system is 
loosely bound. 

PACS numbers: 24.10. Ht; 24.10.Eq; 25.55.Hp 



I. INTRODUCTION 

There is growing interest in using (d,p) reactions as a 
tool to extract single particle properties of nuclei away 
from stability. Since the early studies on 11 Be [![, an in- 
creasing number of transfer experiments have been per- 
formed to obtain detailed structure information. This in- 
cludes experiments at ISOL facilities, performed around 
the Coulomb barrier, such as the recent one on 132 Sn Q, 
but also measurements at higher energy in fragmentation 
facilities, such as the recent systematics on Ar isotopes 
Since the information extracted from these experi- 
ments relies on reaction models, it is critical to validate 
these models and assess their uncertainties. 

There are a number of aspects that need to be consid- 
ered when validating reaction theories. At present, re- 
action theories for (d,p), applicable to nuclei with mass 
A > 10, require a reduction of the many-body problem 
to a few-body problem. A consequence of this first point 
is the introduction of effective interactions between the 
composite constituents, the so called optical potentials 
and single particle binding potential. These are the phys- 
ical inputs to the problem. In addition, there are often 
approximations in solving the few-body problem. It is 
this last aspect that is the focus of the present work. 

For forty years, the tradition has been to use the Dis- 
torted Wave Born Approximation (DWBA) to extract 
spectroscopic information from (d,p) or (p,d) reactions 
[J] . Even though still in use, this method has been chal- 
lenged repeatedly. Just in the last decade, a variety of 
reaction theory studies have been carried out. In [5|, it 
is shown that the choice of the deuteron optical potential 
in DWBA introduces a very large ambiguity which can 
be significantly reduced when using the zero-range ADi- 
abatic Wave Approximation developed by Johnson and 
Soper(ZR-ADWA) [6| built on the nucleon-target optical 
potentials. In ZR-ADWA deuteron breakup is taken into 
account to all orders, while making a zero-range approxi- 



mation for the deuteron. The formalism for a finite range 
version of the adiabatic wave method (FR-ADWA) was 
developed in which is also built on nucleon optical 
potentials. Perhaps even larger than the uncertainties in 
optical potentials, the ambiguity introduced by the single 
particle potentials which describes the many-body over- 
lap function has also been a focus of study @, 0] . Target 
excitation is known to be relevant for specific cases and 
should be considered case by case [1, [l(| . 

To add to the long list of reaction theories, one should 
still consider the continuum discretized coupled chan- 
nel method (CDCC) [l3[ and the Faddeev-type method 
for transition operators, usually referred to as the Alt- 
Grassberger-Sandhas (AGS) method [lj|. Although the 
CDCC method is probably most widely used to describe 
breakup (e.g. '15]), the CDCC wavefunction has also 
been used in the context of (d,p) reactions (e.g. (l6lll7j). 
The AGS method including the Coulomb interaction, 
originally developed for describing few-nucleon reactions 
[18| , has recently been extended to handle (d,p) and (p,d) 
nuclear reactions [r|. The CDCC method takes breakup 
effects into account to all orders without making any fur- 
ther approximations (in some ways, ADWA is an approxi- 
mate version of CDCC) , but Faddeev goes a step further 
in that transfer channels are also included explicitly in 
the expansion. The Faddeev calculations represent thus 
the exact solution to the full three body problem when 
rearrangement channels are present. While both CDCC 
and AGS make less approximations compared to DWBA 
and ADWA, they are computationally demanding. More- 
over, their technical implementations are limited, for ex- 
ample, the present treatment of Coulomb in the AGS 
method has been successfully applied only to nuclei with 
charge Z < 30 [2(3] ■ Finally, and most importantly, these 
two methods are based on a more complex expansion of 
the wavefunction, in such a way that the final cross sec- 
tion depends not only on the many-body overlap of inter- 
est, but on non-trivial interferences between many over- 
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lap functions including overlaps with states in the con- 
tinuum. As a consequence, simpler approximate methods 
are often preferred by experimentalists. 

Recently, a systematic study [llj using the finite 
range version of the Adiabatic wave approximation (FR- 
ADWA) 0, HI] found large finite range effects in (d,p) 
reactions, specially at the higher energies. The trans- 
fer cross section in FR-ADWA depends on the overlap 
function of interest, and not on a complicated superposi- 
tion of terms with many overlap functions. Also impor- 
tant, the method is not computationally expensive and 
is of practical use for the non expert. Nevertheless, one 
needs to appreciate that FR-ADWA relies on a Sturmian 
expansion which is truncated to first order (first order 
should not be confused with DWBA, since the adiabatic 
method relies on an entirely different expansion). This 
first order truncation of the exact solution of the three- 
body problem has never been systematically tested be- 
fore. FR-ADWA calculations in |2l| are used to analyze 
34,36,46 Ar(p,d) data at 33 MeV/u and compared to Fad- 
deev to determine the error in the treatment of the reac- 
tion dynamics. Discrepancies between FR-ADWA were 
found to vary considerably (6-19%). The work in [2l[ 
calls for a better understanding of the range of validity 
of FR-ADWA. 

The aim of the present work is exactly to determine 
the range of validity of ADWA (we drop FR from now 
on, since all ADWA calculations presented here are finite 
range). In order to do that, a systematic comparison 
between ADWA and Faddeev is performed. We cover 
a wide range of beam energies, light and intermediate 
mass nuclei, well bound and loosely bound systems. We 
choose energies for which there is data to ensure that 
the reaction theory performs sensibly, although our fo- 
cus is on comparing two theories starting from the same 
three-body Hamiltonian, and not dwell on detailed com- 
parisons with data. In section II, we provide some key 
aspects of ADWA and Faddeev, in section III results are 
presented and discussed. Finally in Section IV conclu- 
sions are drawn. 



II. THEORETICAL DESCRIPTION 

The starting point for both ADWA and Faddeev is the 
three-body Hamiltonian for n + p + A: 

n 3B = T r + T R + U nA + U pA + V np . (1) 

The interactions between the nucleons and the compos- 
ite target A (U pA and U nA ) should contain an imagi- 
nary term representing the absorption or transfer to other 
channels not explicitly included. We will focus the dis- 
cussion on the reaction A(d,p)B to the ground state (g.s.) 
of B = n + A (note that the formalism for (p,d) in prior 
form is identical to that for (d,p) in post form presented 
here [H). 

Following the work on the zero-range ADWA [6] , John- 
son and Tandy [7] introduced a full finite range ver- 



sion considering the square-integrable Sturmian expan- 
sion. The Sturmian basis is complete within the range of 
the interaction, and is defined through: 

(T r + OiVnJSitf = -EdSitf (2) 

for each given state corresponding to the number of nodes 
i = 0,1,2,.... The inner product between Sturmian 
states is defined by (Si\V np \Sj) = — and at large dis- 
tances, all basis states decay exponentially according to 
the deuteron binding energy. The three-body wavefunc- 
tion is now expanded in terms of Sf. 

oo 

^(r,R) = ^2S i (f) X i(R)- (3) 

i=0 

This form for the three-body wavefunction is then intro- 
duced in the three-body Schrodinger equation for scatter- 
ing H%b^ — Ety which, when imposing the appropriate 
boundary conditions, allows one to calculate Xj(i2). This 
procedure leads to a non-trivial coupled channel equation 
[3 [ll| which could be solved exactly [22j ■ 

Simplicity is recovered when considering only the first 
term in the expansion Eq.Q. 

= S (?)Xo D (R)- (4) 

Then the coupled-channel equation reduces to an optical 
model type equation with distorting potential 

U AD (R) = -(S (r)\V np (U nA + U pA )\S (r)), (5) 

where, apart from a normalization factor, So is the 
ground state wavefunction of the deuteron. 

The three-body wavefunction Eq. Q is now inserted in 
the post-form T-matrix for the (d,p) process [13]: 

T = (4>nAX pB ] \V np + A rem \¥+l) , (6) 

where \ p b IS ^ ne outgoing proton wave, distorted by 
U* B , <f> nA is the wavefunction of the final bound state, 
generated by a binding potential V nAl and the remnant 
term is A rem = U pA - U pB - 

Even if we were not truncating the Sturmian basis, one 
would expect the method to be at its best when remnant 
contributions are small, because then the resulting cross 
section is only sensitive to short distances between the 
neutron and proton, where the basis is complete. How- 
ever, it is not clear that this condition is sufficient when 
the additional simplification Eq.(j4|) is introduced. 



As opposed to ADWA, in the Faddeev approach, elastic 
scattering, transfer, and breakup channels arc treated on 
equal footing. Therefore, it represents the exact solution 
to the problem once a three-body Hamiltonian Eq.((TJ) 
is defined. In the Faddeev method, the wavefunction is 
expanded in an over-complete basis, involving all three 
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Jacobi coordinates [J]. In the AGS method [14j], one 
starts from the Faddcev formalism and arrives at coupled 
integral equations for the transition operators 
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Tpa = (1 - S/3a)G 1 + — 5p~f)tyGoT ya (7) 



7=1 



whose on-shell matrix elements (V^l^al^a) are scatter- 
ing amplitudes and therefore lead directly to the observ- 
ables. In Eq.(0 G = (E + iO - H )- x is the free resol- 
vent, E being the available three-particle energy in the 
center of mass (cm.) system and Hq the free Hamilto- 
nian. The two-particle transition matrix is a solution of 
the Lippmann-Schwinger equation 



ty = Vy + VyGoty 



(8) 



where t> 7 is the potential for the pair 7; we use the odd- 
man-out notation. The channel states \ip~) for 7 = 1,2,3 
are the eigenstates of the corresponding channel Hamil- 
tonian H~ f = Hq + Vj with the energy eigenvalue E; thus, 
I tpy) is a product of the bound state wave function for 
pair 7 and a plane wave with fixed on-shell momentum 
corresponding to the relative motion of particle 7 and 
pair 7 in the initial or final state. 

The Coulomb interaction is an additional complication 
in the AGS method; nevertheless, when only two particles 
are charged and charge numbers Z are not too large, the 
Coulomb interaction has been successfully included [l9l . 
I20I using the method of screening and renormalization 

We solve the AGS equations in the momentum-space 
partial- wave framework [25| . The Coulomb screening ra- 
dius needed for the convergence of the results increases 
with Z thereby increasing also the number of partial 
waves that have to be included. Such calculations are 
not only time consuming but for very high angular mo- 
menta the partial-wave expansion may even become un- 
stable. This limits the application of the technique with 
respect to Z; so far no calculations have been performed 
beyond Z = 30 [2(| ■ Another limitation of the method 
comes from the Pade summation technique for solving 
the AGS integral equation iteratively [26| ; usually it is 
hard to achieve convergence for heavier nuclei at low en- 
ergies. For more details on the numerical techniques we 
refer the reader to ^[H-il. 
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III. RESULTS 

We perform Faddeev AGS and finite range ADWA cal- 
culations for a number of cases that have been repeat- 
edly studied: n Be(p,d) 10 Be(g.s); 12 C(d,p) 13 C(g.s) and 
48 Ca(d,p) 49 Ca(g.s.). Since we are interested in establish- 
ing the range of validity of ADWA, we explore these at 
several beam energies. For nearly all cases we choose 
beam energies for which data is available and we make 
sure that the present theories are able to provide a fair 



FIG. 1: (Color online) Angular distributions for 
11 Be(p,d) 10 Be: (a) E p = 5 MeV, (b) E v = 10 MeV, 
and (c) E p = 35 MeV. 

description of the angular distributions. However, since 
in this work we are concerned with testing ADWA against 
the exact Faddeev AGS theory, we do not include data 
in the plots. 

One of the main inputs to these calculations are the 
pair interactions. In both calculations, a realistic NN 
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interaction is used. Since the ADWA is performed in co- 
ordinate space, and a local form of the NN interaction 
with LS coupling is needed, the Reid interaction [29f is 
used. The more recent momentum based CDBonn [30| 
is straightforward to use in AGS calculations, contrary 
to potentials in coordinate space. These two NN inter- 
actions reproduce the low-energy NN phase shifts, the 
deuteron binding energy, and the quadrupole moment. 
In fact, apart from the short distances, the deuteron den- 
sities produced with Reid and CDBonn are identical, this 
being true also for the 3N and 4N systems. We have ver- 
ified that the use of different NN interactions does not 
induce differences in the predicted (d,p) calculations by 
repeating both FADD and ADWA calculations with a 
simple Gaussian interaction, that reproduces low-energy 
phase shifts and the deuteron binding energy. We con- 
cluded that the transfer cross sections are not sensitive 
to the details of this interaction but the consistency be- 
tween the V np used to obtain the deuteron wavefunction 
and the interaction in the transfer operator (EqJ5]). This 
same conclusion was found in previous works 19]. 

For the nucleon optical potentials we use the global 
parameterization CH89 31]. As all global optical poten- 
tials, the parameters of CH89 are functions of the beam 
energy. In ADWA, when determining ^ ad in Eq.©, half 
the deuteron energy is used for U n A and U p a , whereas 
the exit proton energy is used for the auxiliary potential 
U p b (for calculating XpB and the remnant contribution). 
The same U p a used for ^ ad is used in the transfer opera- 
tor. A real interaction V n A is used to calculate the bound 
state 4> n A- it consists of a Woods-Saxon plus spin orbit 
form with the Woods-Saxon depth adjusted to reproduce 
the experimental binding energy of the nucleus under 
study. In all cases we fix the geometry of the Woods- 
Saxon, the radius r = 1.25 fm and diffuseness a — 0.65 
fm, as well as the spin orbit parameters V so = 6 MeV, 
r so = 1.25 fm, a so = 0.65 fm. Thus, strictly speaking, 
since the neutron interaction is different in the incoming 
and outgoing channel, two different three-body Hamilto- 
nians are used in ADWA. 

In the Faddeev calculation, energy dependence in the 
interactions can introduce orthogonality errors. It is thus 
important to make a wise choice for the energy at which 
U p a is to be calculated. Our standard choice is to fix U p a 
at the energy in the proton channel (FADD), however 
we will also show results using half the deuteron energy 
(FADD2). Concerning the neutron-target interaction, for 
the partial wave where a bound state exists, the same V n A 
interaction is used as for ADWA. For all other partial 
waves, the neutron interaction is U n A, taken from CH89 
at half the deuteron energy [31| . In our standard Faddeev 
calculation (FADD), there are two subtle differences in 
the interactions as compared to ADWA: i) U p a for the 
deuteron channel is obtained at the proton energy, and ii) 
for the partial waves where a neutron bound state exists, 
the neutron scattering potential U n A has no absorption. 
Because we have seen a larger dependence on the transfer 
cross sections to U p a than to U u a in previous calculations 
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FIG. 2: (Color online) Angular distributions for 12 C(d,p) 13 C: 
(a) E d = 7.15 MeV, (b) E d = 12 MeV and (c) E d = 56 MeV 



(see e.g. Fig. 10 of [28j), we believe the uncertainty due 
to i) dominates the overall uncertainty. 

In Fig[T] we present the angular distribution for 
n Be(p,d) 10 Be, at E p = 5 MeV, E p = 10 MeV, and 
E D = 35 MeV. The blue lines are the results obtained 
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FIG. 3: (Color online) Angular distributions for 
48 Ca(d,p) 49 Ca: (a) E d = 19.3 MeV, (b) E d = 56 MeV 
and (c) E d = 100 MeV. 



with ADWA and the red curves are those with Faddeev. 
In order to understand the subtle difference in the inter- 
actions included in the Faddeev calculations compared 



to ADWA, we perform two test calculations: i) a Fad- 
deev calculation using the proton optical potential cal- 
culated at half the deuteron energy (red open circles la- 
beled FADD2) and ii) an ADWA calculation where the 
initial proton distorted wave is calculated with an op- 
tical potential determined at half the deuteron energy 
(purple filled squares labeled ADWA2). The comparison 
between FADD and FADD2 provides an estimate of the 
sensitivity to the choice of energy at which the proton- 
target optical potential is evaluated. If the two methods 
were exact, the comparison between FADD2 and ADWA2 
would provide an estimate of the effect of replacing, in 
the calculation of ^ ad, the appropriate neutron optical 
potential by a real bound state potential in the partial 
wave where a bound state exists. Of course ADWA2 is 
not exact, and therefore the comparison of FADD2 and 
ADWA2 contains both the difference due to the neutron 
optical potential (which should be small) but also a gen- 
uine difference between the treatment of the three-body 
dynamics. Note that if ^ ad were exact, the T-matrix of 
Eq.© should provide cross sections independent of the 
auxiliary potential U p b (used in \pB and A rem ). This 
was recently demonstrated within CDCC [ljj ■ There- 
fore, differences between ADWA and ADWA2 indicate 
inaccuracies in ADWA. 

First of all, there are differences between Faddeev and 
ADWA cross sections at the lower and higher energy. 
At E p = 10 MeV the agreement is perfect. Compar- 
ing FADD with FADD2, we find that only at the higher 
energy is there a dependence in the choice of energy at 
which the proton interaction is calculated, and then Fad- 
deev calculations using the proton interaction at half the 
deuteron energy become slightly closer to the correspond- 
ing ADWA results (compare FADD2 and ADWA2). At 
E p = 35 MeV, ADWA and ADWA2 differ by 11% at 
the peak, which indicates that the truncation to the first 
term in ADWA is insufficient at these energies. 

In Fig[5]the same observable is shown for 12 C(d,p) 13 C 
at E d = 7.15 MeV, E d = 12 MeV and E d = 56 MeV. 
Differences between the Faddeev and ADWA are small- 
est at E d = 12 MeV. As for u Be, only at the highest 
beam energy are the results sensitive to the choice of en- 
ergy at which the proton optical potential is defined. In 
this case, calculating the final proton distorted wave with 
UpB evaluated at half the deuteron energy slightly im- 
proves the agreement (compare ADWA2 versus FADD2). 
This is not to say that adiabatic calculations for (d,p) 
should be performed using U p b at half the deuteron en- 
ergy. The standard ADWA makes the correct choices. 
Here ADWA2 only helps to understand which part of the 
disagreement with Faddeev comes from the slight differ- 
ences in the proton interaction. 

Our last case 48 Ca(d,p) 49 Ca is shown in Figf51 An- 
gular distributions were calculated at E d — 19.3 MeV, 
E d = 56 MeV and E d = 100 MeV. The comparison be- 
tween ADWA and Faddeev follows the same trend as in 
the previous cases: good agreement at around 10 MeV/u, 
with large deterioration at very high energies. It was not 
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possible to obtain Faddeev solutions for energies lower 
that 19.3 MeV due to lack of convergence of the Pade 
summation technique. Oddly, for the reactions at 56 
MeV, while ADWA showed a very strong dependence 
on the choice of U p b, the Faddeev results only exhibit 
a dependence on the proton optical potential around the 
second peak. At 100 MeV, the angular distribution is ex- 
tremely sensitive to the choice of the proton energy used 
to determine the proton optical potential parameters for 
both ADWA and FADD. This introduces large ambiguity 
in the comparison of ADWA and Faddeev, at 100 MeV. 

Overall, we find the differences between ADWA and 
Faddeev to be below 10 % except for the lowest and high- 
est energies. Our results suggest that there is an optimum 
beam energy (around 10-20 MeV/u) where ADWA is at 
its best, differing from Faddeev by only a few percent. 
Apart from the 100 MeV 48 Ca(d,p), changes introduced 
in the shape of the angular distributions by the adiabatic 
approximation are small. A quantitative summary of our 
results, plus additional detail in our calculations, is given 
in Table |TJ Here Aj?_,4d is the percentage difference be- 
tween cross sections from FADD and ADWA, relative to 
the FADD cross section determined at the first peak of 
the distribution (corresponding to 6). To study the effect 
of the choice of the energy at which the proton interac- 
tion is determined, we also show Ap-Fi, the percentage 
difference between Faddeev using the energy in the pro- 
ton channel to determine U p a (FADD) and that using 
half the deuteron energy (FADD2). This represents the 
ambiguity in the three-body Hamiltonian for the Fad- 
deev calculations, and can be used as an uncertainty in 
the comparison. It is not important for low and interme- 
diate beam energies, and therefore in this beam-energy 
region, the agreement between Faddeev and ADWA is ro- 
bust. For completeness, the identical quantity for ADWA 
Aad-ad2 is also shown (where the standard ADWA is 
compared with that where the proton distorted wave 
is calculated with CH89 using half the deuteron energy 
ADWA2). As the exact T-matrix should be independent 
on the choice of the auxiliary potential, Aad-ad2 pro- 
vides an internal warning sign that the approximation in 
ADWA is not adequate. Consistent with the compari- 
son with Faddeev, the effect becomes significant for the 
higher beam energies. 

In Section II, we pointed out that the Sturmian expan- 
sion is expected to be good in the region of V np but not 
for large neutron-proton distances. For this reason, we do 
not expect ADWA to work when the contribution of the 
remnant term in Eq.([S]) is large. The last column of Table 
II contains the percentage difference between the original 
ADWA and the ADWA calculation neglecting the rem- 
nant term, relative to the original ADWA. The only case 
for which the remnant contribution is large is for the high 
energy reaction with 11 Be and then, indeed, ADWA per- 
forms poorly in comparison with Faddeev. However in 
the other cases, it is not the remnant contribution that 
is responsible for the disagreement ADWA versus Fad- 
deev. Since the largest disagreement comes at the highest 



TABLE I: Percentage differences between the differential cross 
section at the first peak of the distribution for the various 
formulations: Af-ad comparing Faddeev with Adiabatic, 
Af-up the effect of changing the energy at which the proton- 
target interaction is calculated in the Faddeev, and finally 
AAD~rem the effect of the remnant term in the adiabatic. 
Also given is the beam energy E (in MeV) and the angle 9 at 
which the percentage difference was calculated (in degrees). 



reaction E 9 


Af-AD Af-F2 AaD-AD2 &-AD-rem 


11 Be(p,d) 5 1 
10 1 
35 1 


-22.90 0.17 -0.10 1.36 
-1.12 0.42 0.44 -1.13 
18.50 15.25 10.9 26.82 


12 C(d,p ) 7 20 
12 15 
56 1 


-9.08 1.61 -0.42 -3.74 
-2.74 2.18 -0.50 -3.92 
-14.26 13.95 17.5 3.34 


48 Ca(d,p) 19 8 
56 1 
100 


-6.12 7.16 -0.64 -0.74 
-3.05 0.76 46.4 0.57 
51.0 72.9 -42.7 9.73 



beam energies, this suggests that the main source for the 
ADWA disagreement with Faddeev comes from the trun- 
cation of the Sturmian expansion Eq.([5]). Eq.Q implies 
that all xt D {R) behave as the elastic Xo D (R)- However, 
at higher beam energies, large excitation energies in the 
deuteron channel are expected to affect the shape of the 
adiabatic distorted waves Xi(R)- 

An attempt to systematize the results is done in FigfJ] 
where the ratio of Faddeev cross sections at the first 
peak of the angular distribution (9 from Table UJ over 
the corresponding ADWA cross sections is plotted as a 
function of EJ™ /V c , the cm. energy of the deuteron 
over a simple estimate of the Coulomb barrier 
V c = 1.44Z/(1.2A 1 / 3 ). We introduce an error bar cor- 
responding to the ambiguity in the choice of the energy 
at which the proton interaction is calculated in the Fad- 
deev (2nd column in Table fl]). Figure 2] provides a good 
illustration of the validity of ADWA for intermediate en- 
ergies. For the lowest energies, differences increase up to 
20% and they become particularly large for the highest 
energies where the Hamiltonian ambiguities are also the 
largest. 

In addition there should also be an error due to the fact 
that in the Faddeev approach, the n-target interaction in 
the deuteron channel has no absorption for the partial 
wave in which there exists a bound state. Based on pre- 
vious experience, we expect this effect to be smaller than 
the effect on the choice of U p a [Hj]. Additional Faddeev 
calculations for the lighter targets, introducing energy de- 
pendence in the n-target interaction in the partial wave 
corresponding to the bound state, suggest the ambiguity 
on U n A is only comparable to Ap~F2 for 11 Be(p,d) at 35 
MeV. 

We have also explored the dependence on the angu- 
lar momentum I of the neutron bound state as well as 
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FIG. 4: (Color online) Ratio of Faddeev prediction for the 
cross section at the first peak of the angular distribution ver- 
sus the Adiabatic counterpart plotted in term of the deuteron 
energy in the cm. over the Coulomb Barrier. 



its binding energy — S n . For this purpose, we repeat 
the calculations for our best case 11 Be(p,d) 10 Be at 10 
MeV. The neutron bound state in 11 Be has I = and 
S n = 0.5 MeV. We first repeat the calculations chang- 
ing the angular momentum to I — 1,2,3 by adjusting 
the depth of the neutron bound state interaction V u a 
to reproduce the binding S„ =0.5 MeV, while keeping 
all other interactions fixed. We find that as I increases, 
the difference between FADD and ADWA also increases: 
A(F - AD) = 1% for I = 1; A(F - AD) = 15% for I = 2 
and A(F - AD) = 30% for I = 3. 



Next we repeat the calculations fixing I = and chang- 
ing the depth of V n A to reproduce S„ = 2.0,4.0 and 6.0 
MeV. Again, all other interactions are unchanged. The 
percentage difference relative to the Faddeev result is: 
A(F - AD) = 1% for S n = 0.5 MeV, A(F - AD) = 13% 
for S n = 2.0 MeV, A(F - AD) = 18% for S„ = 4.0 
MeV and A(F - AD) = 23% for S n = 6 MeV. In these 
calculations, remnant effects do not change significantly. 
However, the deeper the neutron initial state, the more 
delocalized is the n — p system in the continuum, and 
one might expect the first term of the Sturmian expan- 
sion Xq D {R) to have the wrong radial behavior. Whereas 
with increasing / the ADWA diverges from the Faddeev 
in a rapid manner, the differences with increasing binding 
appear to reach a saturation point around 20 %. These 
results explain the findings in [2lJ. In [2lJ, the discrep- 
ancy between ADWA and Faddeev was found to be much 
larger for 36 Ar(p,d) than for 34 Ar(p,d) and 46 Ar(p,d). 
The hole state created by pulling out a neutron from 
36 Ar has both, large I as well as large S n . 



IV. SUMMARY AND CONCLUSIONS 

The finite range adiabatic model provides a practical 
method to analyze (d,p) and (p,d) reactions that goes 
well beyond DWBA. One of the great advantages is that 
it uses nucleon optical potentials, rather than deuteron 
optical potentials which are far more ambiguous. How- 
ever it is by no means the exact solution and in this 
work we quantify the errors introduced by comparing it 
to exact Faddeev AGS calculations. We choose three em- 
blematic reactions n Be(p,d) 10 Be(g.s.), 12 C(d,p) 13 C(g.s.) 
and 48 Ca(d,p) 49 Ca(g.s.) and span a large range of beam 
energies. Overall, we find the agreement for reactions 
with beam energies around 10 MeV/u to be better than 
10%. These results bear important implications for the 
reactions being measured at Isotope-Separator-on-Line 
facilities as well as for the science program at the future 
Facility for Rare Isotope Beams. The deviation of ADWA 
from Faddeev increases with the angular momentum of 
the neutron bound state as well as with the separation 
energy of the neutron. The dependence on I was found 
to be stronger than that on the separation energy. 

The comparison with Faddeev is limited in several 
ways. On one hand there are limitations with fixing the 
three-body Hamiltonian. On the other hand there are 
technical issues that limit the number of cases that can 
be studied. To address the first point, we have explored in 
detail the effect of different choices for the proton target 
interaction, namely on the proton energy used to deter- 
mine the optical potential. The differences obtained are 
then used as a systematic uncertainty in the comparison. 
They are only important at the higher energy, exactly 
where the performance of ADWA is at its worse. In this 
sense, the conclusion that ADWA fails for energies larger 
than E c d m /V c w 12 is not robust. 

As mentioned in Section II, there are several techni- 
cal limitations in the present implementation of the Fad- 
deev equations, namely the AGS method in momentum 
space. In this respect we found 48 Ca(d,p) at 19.3 MeV 
to be most challenging and close to the limit of present 
capabilities in obtaining converged results. 

Since the finite range adiabatic model for (d,p)/(p,d) 
at intermediate energies appears to be well suited to de- 
scribe the few-body reaction dynamics for many cases of 
interest, is easy to use and is not computationally inten- 
sive, it is strongly desirable that the standard reaction 
codes incorporate this option in a user-friendly manner. 

In this work we did not consider target excitation, how- 
ever it is understood that there will be many applications 
in which target excitation is an integral part of the reac- 
tion mechanism. It would thus be worthwhile to extend 
ADWA to include the inelastic channels in the formula- 
tion consistently. 
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